Persistence of equilibrium states in an oscillating double-well potential 
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We investigate numerically parametrically driven coupled nonlinear Schrodinger equations mod- 
elling the dynamics of coupled wavefields in a periodically oscillating double-well potential. The 
equations describe among other things two coupled periodically-curved optical waveguides with 
Kerr nonlinearity or horizontally shaken Bose-Einstein condensates in a double- well magnetic trap. 
In particular, we study the persistence of equilibrium states of the undriven system due to the 
presence of the parametric drive. Using numerical continuations of periodic orbits and calculating 
the corresponding Floquet multipliers, we find that the drive can (de)stabilize a continuation of an 
equilibrium state indicated by the change of the (in) stability of the orbit. Hence, we show that 
parametric drives can provide a powerful control to nonlinear (optical or matter wave) field tunnel- 
ing. Analytical approximations based on an averaging method are presented. Using perturbation 
theory the influence of the drive on the symmetry breaking bifurcation point is discussed. 



I. INTRODUCTION 

Parametric drives, i.e. external drives that are periodic 
in time and depend on the system variables, have been 
used as a means to control and maintain a system out- 
of-equilibrium [I] |2] . The generation of standing waves 
when a liquid layer is subjected to vertical vibration, 
known as Faraday waves, is among the classical stud- 
ies of parametrically driven instabilities [3 . The vertical 
vibration can sustain spatially localized, temporally oscil- 
lating structures, commonly referred to as oscillons, such 
as in granular materials [4, 5 , Newtonian [6 , 7 and non- 
Newtonian fluids [8 . Relatively recently, ac parametric 
drives have been predicted and shown to be able to sus- 
tain localized waves in a linear Schrodinger system. In 
the undriven case the excitations would simply disperse 
[9]. This effect is referred to as dynamic localizations. 
Recently, it has been suggested theoretically [TOl |TT] and 
shown experimentally [TU |13] that periodically curved 
optical waveguide arrays can be an ideal system for real- 
izations of dynamic localizations. Such a localization has 
been used in quantum physics in the context of Bose- 
Einstein condensates to reduce or even completely sup- 
press quantum tunneling of particles trapped in a po- 
tential well by shaking the potential back and forth (see 
the review [TT). The tunneling suppression in a series of 
potential wells using the method has been shown experi- 
mentally p!5HT7] . The method has been proposed as well 
as a powerful tool to manage the dispersion of, e.g., an 
atomic wavepacket [18] and nonlinear gap-solitons [19]. 

Complementing the seminal finding of dynamic local- 
ization, coherent destruction of tunneling predicted in 
[2Q[ [2T] has become an important phenomenon in the 
study of quantum dynamical control. While dynamic 
localization occurs in infinite dimensional systems, sys- 
tem boundaries play an important role in the coherent 
destruction of tunneling. The internal relationship be- 
tween the two phenomena was discussed recently in [22 
where it was shown that both phenomena can be inter- 
preted as a result of destructive interference in repeated 



Landau-Zener level-crossings. A recent experiment on 
strongly driven Bose-Einstein condensates in double- well 
potentials formed by optical lattices provided a direct ob- 
servation of dispersion suppression in matter waves [23]. 
In the context of optical physics, the idea of tunneling 
destruction in couplers, i.e. dual-core waveguide arrays, 
was proposed in [24 and later implemented experimen- 
tally in [25l[26]. A sketch of the physical setups is shown 
in Fig. [Tj 
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FIG. 1: Sketch of (a) periodically-curved couplers, (b) Bose- 
Einstein condensates in a periodically oscillating double-well 
potential. 

It is important to note that the above phenomena are 
in the linear regime (see a recent review on theoretical 
and experimental advances in modulated photonic lat- 
tices that provide controls on fundamental characteris- 
tics of light propagations [27 ). The parameter values for 
tunneling suppression are isolated degeneracy points of 
the quasienergies. When nonlinearity is present in the 
system, the tunneling can be suppressed in a relatively 
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wide interval of parameter values [28]. It is also inter- 
esting that periodic driving may enhance tunneling as 
opposed to preventing it [29l [30] , which has been shown 
experimentally recently [31]. 

Despite many important works mentioned above, it 
is necessary to note that so far there is no study on 
the effects of parametric drives from a dynamical sys- 
tem point of view. In that regard, the present work 
aims at studying the effects of nonlinearity on the afore- 
mentioned tunneling enhancements or suppression. In 
particular, we consider the persistence, i.e. the existence 
and stability, of equilibrium states when the potential 
is periodically oscillating. As discussed in the theoreti- 
cal [32, 33 and experimental [34-36, studies of Joseph- 
son tunneling of Bose-Einstein condensates in a double- 
well potential, the dynamics of the macroscopic quantum 
density and phase difference between two condensates is 
described by two coupled nonlinear ordinary differential 
equations. Depending on a control parameter, one can 
obtain qualitatively different dynamical behaviors that 
can be explained well by classical bifurcation theory. In 
general, there are three types of fixed-points in the phase- 
portraits. Using notations in the limit of no coupling be- 
tween the condensates, the equilibria can be symbolically 
written as [+,+], [+,0], and [+, -] where and "-" 
represent and 7r-phase of the wavefunctions in the two 
potential minima, respectively, and "0" the case when 
the density vanishes in the well. The second equilibrium 
state is also commonly referred to as a self-trapped state 
because the density difference between the wells oscil- 
lates about a nonzero value. The effects of parametric 
drives, i.e. an oscillating potential, on such equilibria are 
studied numerically in the present work. It will be shown 
that the parametric drive in general does not influence 
the existence of an equilibrium, which becomes a rela- 
tive periodic orbit due to the drive, but indeed may alter 
the stability. In particular, we will show that the self- 
trapped states, i.e. no-tunneling states, can be broken by 
the drive creating, possibly chaotic, tunneling between 
the two wells, i.e. a destruction of self-trapping with the 
relative phase of the fields oscillating decoherently. In 
addition to that, we will also show that the symmetric 
[+, +] -state which is unstable beyond a critical norm can 
be made stable by the drive, which we call a coherent con- 
struction of tunneling. The antisymmetric [+, — ]-state, 
which is generally robust, can also be destabilized by the 
parametric drive. 

The paper is presented as follows. In Section II, we 
discuss the governing equation and the numerical meth- 
ods used in the paper. In the following section, the de- 
struction of self-trapped states due to a parametric drive 
is presented. We show that the destruction is caused 
by a period-doubling bifurcation. In Section IV, we 
study the persistence of symmetric states. We show that 
parametric drives are able to stabilize or destabilize the 
states. A similar numerical discussion for the antisym- 
metric state is presented in Section V. The state which 
is always stable in the undriven case can become unsta- 



ble due to parametric drives. Analytical calculations are 
presented in Section VI. We discuss the (de) stabilization 
intervals observed in the previous sections using an aver- 
aging method. The symmetry breaking bifurcation point 
where the asymmetric and symmetric states merge is 
shown using perturbation theory to be affected by para- 
metric drives. Conclusions are presented in the last sec- 
tion. 



II. GOVERNING EQUATION AND 
NUMERICAL METHODS 

Using a tight-binding approximation, the parametri- 
cally driven wavefunctions due to an oscillating potential, 
i.e. in periodically oscillating waveguides in nonlinear op- 
tics or horizontally shaken double-well magnetic traps in 
matter waves, are described by (see, e.g., [TTJ |T3l |T8] ) 

iui = 5 \uif Ui — qui -\- ce~*^°^^^ii2, 
iu2 = S \u2f U2 — qu2 -\- ce'^^^^^^ui. 

In the context of nonlinear optics, Uj is the optical field 
in the jth waveguide, c > is the waveguide coupling 
coefficient (in units of 1/mm), the dot is a derivative 
with respect to the propagation direction z, q is the light 
propagation constant and (5 > is the nonlinearity co- 
efficient (l/(W-mm)). The defocusing case ^ < can 
be obtained immediately due to the staggering trans- 
formation Uj {—lyuj. The parametric drive is rep- 
resented by the function Xo{z) = {nsO:/h)xo{z), where 
a is the separation distance between the waveguides in 
units of /im, rig is the substrate refractive index, xo{z) 
describes the physical periodic curving profile and h is 
the inverse of the light wavenumber. In this work, we 
take xo{z) = auj sm{ujz) with amplitude a in units of mm 
and waveguide curvature wavenumber u. The constant, 
time-independent quantity = \ui\'^ + |ii2p associated 
with the field power is referred to as a norm herein. The 
parameters are taken in the vicinities of those used in 
[28] . The onset of, e.g., the existence and the stability 
of a solution will certainly depend on the parameter val- 
ues, nevertheless the results presented herein are qualita- 
tively generic. Without loss of generality, one can scale 
6 = UJ = 1. When a = 0, Eq. ^ is integrable. 

To study the persistence of an equilibrium in the pres- 
ence of a parametric drive, we solve the governing equa- 
tion ([T]) for periodic orbits. We therefore seek solutions 
satisfying Un{0) = u^iT)^ where T = 27r is the oscillation 
period. In the absence of a drive, equilibrium solutions of 
the equations clearly fulfil the relation. In the presence 
of drives, we look for the continuations of an equilibrium 
point by using shooting methods in real space or alge- 
braic methods by discretizing the propagation direction 
variable. The latter method is more useful than the first 
when there is a bifurcation in the numerical continuation. 
In that case, we use a pseudo-arclength continuation al- 
gorithm [37H39]. 
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When a periodic orbit, say Un{z), is obtained, we also 
examine its stability by calculating its Floquet multi- 
pliers, which are eigenvalues of the monodromy matrix. 
This is obtained from solving a linearized equation about 
the solution Un{t) 

Oil =5{2\Ui\^ui^Ulu\ \ -g^ii +ce-^^o(^)iX2, 
m2 = 5(2\U2\^U2^Uiul^ -qu2^ce'^'^'^ui. 

The linear system is integrated using a Runge-Kutta 
method of order four or a symplectic method. The 
ith column of the monodromy matrix M is the vec- 
tor [Re(Mi(T)),Re(M2(T)),Im(Mi(T)),Im(M2(T))]T, 
that corresponds to the initial condition 
[Re{ui (0)), Re(ii2(0)), lm(iii(0)), Im(^i2(0))]^ that is 
equal to the ith column vector of the identity matrix I4. 
A solution is stable when all the multipliers are on the 
unit circle. Note that here it is not necessary to have 
a pair of eigenvalues at (+1). The (in) stability result 
obtained from calculating the monodromy matrix is also 
confirmed by integrating Eq. ([T]). 

Writing Uj = \uj\e'^^^^^ it is convenient to represent 
a solution with its population imbalance A and phase- 
difference between the light fields described as 

A = (|Mi|2-|u2|2)/iV, e = cP2-cPi. (3) 

One can show that A and satisfy parametrically driven 
sine- Gordon equations 

A = 2 cVl — A^ (sin xq cos 6 — cos xq sin 6>) , 
6 = (cos xo cos e + sin xq sin (9) - ANS. 

When there is no drive, i.e. xq = 0, we obtain the equa- 
tions derived in [32 . In that case, Q will have at most 
three fixed points given by (^, A) = (0,0), (±7r,0), and 

(0, ±^1 — (2c/((5A^))^), which correspond respectively to 
the symmetric, antisymmetric, and asymmetric state. It 
is clear that the asymmetric state pair only exists when 
c < SN/2. The symmetric and antisymmetric states in 
the focusing case correspond respectively to the antisym- 
metric and symmetric states of the defocusing case due 
to the staggering transformation. 

III. DESTABILIZATION OF ASYMMETRIC 
SELF-TRAPPED STATES 

First, we consider the effect of a parametric drive on 
the asymmetric ground state. Let us set the propaga- 
tion constant, e.g., q = 2. When a = 0, the bifurcation 
diagram of the stationary [+, +] state is shown in Fig. 

The symmetric [+, +] -state corresponds to branch BD. 
For a sufficiently small norm, the state is known to be sta- 
ble (see, e.g., [40l[4T] and references therein). When one 
decreases the parameter c further, there will be a critical 



E 




^"x A 




0.7 0.8 


0.9 


1 1.1 




C 






(a) 






(b) 



FIG. 2: (a) Bifurcation diagram of equilibrium states in 
straight couplers a = 0. Shown is the norm of the states 
against the coupling parameter c with q = 2. Dashed line 
indicates unstable solutions, (b) Phase-portraits of ([T]) with 
a = 0, = 2 and c = 10/11. 



norm above which the symmetric state is no longer stable. 
The equilibrium loses stability at a pitchfork (symmetry 
breaking) bifurcation with an asymmetric state leading 
to " macroscopically quantum self-trapping" in the con- 
text of matter waves [32l|33], i.e. point A in the figure. 
The state corresponding to branch AE can be viewed 
as the [+,0]-state. In Fig. [2|b), we show the phase- 
portraits of the coupled equations ([T]) in the (6>, A)-plane 
with N = 2 and c = 10/11. One can see that the sym- 
metric state represented by (0, 0) is indeed unstable while 
there is a pair of stable fixed points lying on the verti- 
cal axis ^ = 0, which is the asymmetric ground state. 
The BC branch in Fig. |2|a) corresponds to a periodic 
orbit encircling both the stable fixed points with = 
in Fig. [2|b). The branch FG corresponds to a periodic 
orbit encircling only one of the stable fixed points. Both 
branches do not correspond to equilibria, but are shown 
for the completeness of the analysis later. 

Next, we consider the presence of a parametric drive 
in the system. Rather than showing the phase portraits 
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FIG. 3: (Color online) Stroboscopic plots of the system ([T]) 
with AT = 2 and (a) c = 10/11 and (b) c = 2/3. Dots with 
the same color are obtained from the same initial condition. 



in the A)-plane with a continuous time, it becomes 
more convenient to represent the solution trajectories 
in Poincare maps (stroboscopic plots at every period 
T = 27r). Note that in a Poincare map a periodic or- 
bit will correspond to a fixed point. In the recurrence 
map, a stable periodic orbit will therefore correspond to 
an elliptical fixed point encircled by closed regions (is- 
lands). 

Setting a = 0.1, we show in Fig. |3] the recurrence maps 
of the system for norm N = 2 with two different values of 
c, i.e. c = 10/11 and 2/3, in the (^, A)-plane. The plots 
are obtained from various sets of initial conditions using 
direct numerical integrations of the governing equations 
([T]). It is important to mention that in the first case, there 
is no asymmetric trapped state. The island with ^ = in 
the figure is actually a tunneling state (A(t) is not sign- 
definite), that is locked to and oscillates coherently with 
the drive. It is also clear that there are rather chaotic 
oscillations between the fields in the two potential min- 
ima. Interestingly, by changing the coupling constant, 
one now has self-trapped states shown in Fig. [sf^b) by a 
pair of small islands centered at ^ = 0. The question is: 




FIG. 4: Bifurcation diagram of equilibrium states in periodi- 
cally curved couplers with a = 0.1. The inset shows Floquet 
multipliers of the asymmetric state in the complex plane for 
c= 10/11. 



what happens with the states in the first case? We con- 
jecture that the two topologically different maps in Fig. 
|3]are determined by the persistence of such fixed points. 
To illustrate our conjecture, we have solved the govern- 
ing equations ([T]) for periodic orbits. The corresponding 
bifurcation diagram of Fig. [2|a) when a = 0.1 is shown 
in Fig. [4] where a rich bifurcation structure is evident. 

Comparing Figs. [2|a) and|4j one can see that there is 
also branch merging and splitting. Branches FG and BC 
in Fig. [2|a) split into FjGj and BjCj, j = 1,2 in Fig. 
|4] More importantly, the branch FiGi merges with a 
segment AF that corresponds to an asymmetric ground 
state, which now becomes unstable. The norm used in 
Fig. [3]^ a) belongs to this unstable region, explaining why 
there is no asymmetric trapped state in the Poincare sec- 
tion. In the inset of Fig. |4j we display the Floquet mul- 
tipliers of the state showing that it suffers from an in- 
stability due to a pair of multipliers leaving at —1, i.e. a 
period-doubling bifurcation. The question whether sub- 
sequent period-doubling bifurcations occur in the system 
is left for further future investigations. 

The presence of only one family of islands in Fig. [sja) 
centered at ^ = 0, i.e. without its counterparts in the 
lower half-plane, can also be explained using the bifurca- 
tion diagram in Fig. [4j When a = 0, one of the islands 
corresponds to a periodic orbit along the branch BC. 
As the branch splits into unstable BiCi and stable B2C2 
branches, it is expected that there will be only one family 
of stable coherent periodic orbits as seen in the Poincare 
sections. The splitting itself is expected from the Floquet 
multipliers of the BC branch as they are all at +1 (not 
shown here). This is also the case with branch FG. 

The destabilization of the asymmetric state is rather 
generic, including the case of small N. One difference is 
that we only observed unstable solutions and did not see 
any branch splitting along the asymmetric state branch 
unlike in Fig. [4] 
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FIG. 5: (a) The norm of symmetric states as a function of the 
amphtude drive a for c = 0.1 with q = 2 (upper) and q = I 
(lower). Dotted Une indicates unstable solutions. The insets 
show the Floquet multipliers of two solutions along the upper 
curve at a = 2.5 and a = 3.5. (b) The time dynamics of the 
initial condition ui = U2 with c = 0.1, N = 2 and a = 3.5. 
The inset depicts the trajectories of the dynamics in time in 
the (6*, A)-plane. 



IV. PERSISTENCE OF SYMMETRIC STATES 

We have seen in Section HIl that in the undriven case 
when the solution norm exceeds a critical value, the sym- 
metric states lose stability to asymmetric states, see Figs. 
[2] It is natural to question whether it is possible to stabi- 
lize unstable symmetric states using the parametric drive. 
We show in Fig. [5] that it is indeed the case. 

In Fig. [sja) we depict the norm of symmetric states 
as a function of the drive amplitude a for two values of 
q. As a particular choice for clarity, we consider a sym- 
metric state with a relatively large norm and small cou- 
pling such that when there is no drive the wave fields are 
weakly coupled and rather strongly unstable. It is inter- 
esting to note that upon increasing the parameter a, the 
continuation of the undriven symmetric state becomes 
stable. The stabilization is generic as we also observe 



the same phenomenon for other values of q with the sta- 
bilization threshold almost independent of q as shown in 
the lower curve in Fig. [sja). In the inset of the figure, 
we show the Floquet multipliers of the states at a = 2.5 
and a = 3.5 corresponding to values before and after the 
stabilization threshold. For the second case, the periodic 
orbit is stable shown by the multipliers that are all on 
the unit circle. 

To see whether the stabilization observed in the con- 
tinuation of the periodic orbits above can be viewed from 
direct numerical integrations of the governing equations 
([1]), we show in Fig. [sjb) the time dynamics of the initial 
condition m = U2 dit z = showing a symmetric state 
with coherent relative phase between the fields. Hence, 
we have obtained a coherent construction of tunneling in 
the region when symmetric (in-phase) states are unsta- 
ble. 




z 



(b) 

FIG. 6: (a) The same as Fig. [sja), but for stable symmetric 
states (when undriven) with c — 0.4. The inset shows the 
Floquet multipliers of a solution along the upper curve at 
a = 2.5. (b) Typical dynamics of unstable symmetric states 
with a — 2^ N — 0.7 and 1x1 (0) = 1/^2(0). The inset shows the 
dynamics in the (^, A)-plane. 

In addition to stabilization, we have also studied the 
possibility of using a parametric drive to destabilize sym- 
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metric states that are stable in the undriven case. Our 
study is summarized in Fig. |6|a) for c = 0.4, where we 
show that the continuation of an undriven symmetric 
state becomes unstable in an interval of driving ampli- 
tude. As for the dynamics of the instability, in Fig. ^h) 
we depict a possible manifestation of the destabilization 
for ui{0) = 2x2(0) and N = 0.7. For this value of power, 
the state is stable when undriven as it is below the crit- 
ical norm for a pitchfork bifurcation. When the state is 
driven, the relative phase is no longer localized and the 
tunneling is completely suppressed, which can be called 
a nonlinear coherent destruction of tunneling similar to 
that discussed in [28 . Note that the instability window 
observed in Fig. [6[a) is relatively wide. We believe that 
the finite range for the nonlinear coherent destruction of 
tunneling reported in [28^ is due to an instability of the 
periodic orbits. 



V. DESTABILIZATION OF ANTISYMMETRIC 
STATES 

Finally, we consider the antisymmetric [+, —] -state, 
which is stable in the undriven case. The state is gen- 
erally robust. Nevertheless, in a similar fashion to the 
symmetric state in the previous section, we show in Fig. 
[7| that antisymmetric states can also become unstable 
when driven above a threshold value. The instability do- 
main is finite similar to the symmetric states discussed 
in Section IV. When one increases the drive amplitude 
further, there will be other instability windows. 

In Fig. [7|^b), we show the typical dynamics of unstable 
[-h, —] -state, where we also obtain a nonlinear coherent 
destruction of tunneling similar to Fig. [ojb). 



VI. ANALYTICAL RESULTS: STABILIZATION, 
DESTABILIZATION AND SYMMETRY 
BREAKING BIFURCATION POINTS 

In this section we discuss the effects of parametric 
drives on the continuation of equilibrium states analyti- 
cally. 

A method commonly used in almost all of the previous 
theoretical studies is to average the governing equation 
([1]). The approach is justified when the driving frequency 
is large. In our scaling, that corresponds to small norm 
solutions. The averaged equation is given by 



S\uj\ Uj - quj ^ cJo{a)u2-j, j = 1,2, (5) 



where Jq is the Bessel function of the first kind. Hence, 
we obtain a coupler with an effective coupling constant 
cJo(a). 

As discussed in Section I and II above, there is a critical 
threshold at which the symmetric and asymmetric states 
change stability, i.e. the symmetry breaking bifurcation 
point c = N/2. Note that Jo (a) oscillates about with 
the oscillation amplitude decreasing with a. Therefore, 
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FIG. 7: (a) The same as Fig.|5ja), but for the antisymmetric 
state with c = 0.1. The inset shows the Floquet multipliers 
of a solution along the upper curve at a = 2.5. (b) The same 
as Fig. |6jb) for the antisymmetric mode with N — 2 and 
a 2.5. 



in the context of the average equation the only possible 
way that the continuation of an equilibrium state changes 
stability is when the effective coupling crosses ±A/'/2. 

For the asymmetric state discussed in Section III, it is 
expected that the state becomes unstable in the interval 
of a where cJo{a) > N/2. For the symmetric state stud- 
ied in Section IV, the state would become stable when 
cJo{a) > N/2 or cJo{a) < 0. The latter is due to the stag- 
gering transformation as the state effectively becomes an 
antisymmetric one when cJo(a) changes sign. Finally for 
the asymmetric state analyzed in Section V, the state is 
expected to become unstable when > cJo(a) > —N/2 
due to the staggering transformation. The analytical ex- 
planations of the stability switching above are rather in 
good agreement with the numerical results. 

Despite the agreement, the averaged equation fails in 
describing several phenomena behind the stabilization 
and destabilization of the original governing equation, 
such as the period-doubling bifurcation that causes the 
destabilization of the asymmetric state instead of the 
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state ceasing to exist and the branch spUtting. Another 
failure of the averaged equation is in predicting the in- 
fluence of parametric drives on the symmetry breaking 
(pitchfork) bifurcation point. Because Jo (a) < 1 for 
a > 0, at the bifurcation point of the undriven case, i.e. 
c = the effective couphng cJo{a) is less than half 

of the solution norm, which implies that the symmetric 
state should be unstable when the system is parametri- 
cally driven. Nevertheless, from Fig. [4] we obtained that 
the bifurcation point occurs at a larger norm, which con- 
tradicts the averaged equation. To resolve the issue, we 
used a perturbation expansion. 

We consider the equivalent equation Q. Studying the 
continuation of the symmetric state (A,^) = (0,0) in 
the presence of a parametric drive with small amplitude 
\a\ <C 1 at the bifurcation point of the undriven case 
c = N/2, we take the following expansion 



A = aA^i) + 0{a^), = aO^^^ + ©(a'), 
which upon substitution into Q yields 

0. 



^A^i) = -7V^(i)+7Vsinz, ^6^^^ 
dz dz 



Looking for periodic solutions, we obtain 
A^^) = -Ncosz, e^^^ = 0. 



(6) 



(7) 



Next, we study the stability of the periodic solution 
above. The linearisation of Q about periodic solution 
A and 6 is 



X — 



2cA 



(sin xq cos 6 — cos xq sin 6) x 



2 c^/Y— A^ (sin Xq sin + cos xq cos 0) 

(t^2- + 1^ (cosxo cos^ + sinxo sin^) X 
(— cos Xq sin 6 + sin xq cos 0) y. 

(8) 



-NSx 



2Ac 



We will use the solution of the equation to construct a 
Floquet matrix that will determine whether or not the 
periodic solution is stable. It is then natural to expand 
the variables x{t) and y{t) in series as A and above, i.e. 



X = +a2a;(2) +0(a3), y = y(°^ +a^y(^HO{a^). (9) 



Substituting the expansions into the linearized equation 
([8|, from terms of 0{1) we obtain 



fc(0) = -iVy(0),y(0)=0. 



(10) 



Prom terms of 0{a^) we have 



-Ny^ 



2) 



sm 



y(2) 



^sin(2z)yW. 



+ f {N^ (cos(2^) + 1) + 1 - cos(2z)) 
f (3iV2 - 1) (cos(2^) + 1) a;(0) - ^ si: 

-S ( |To| are subject to the initial condit 
either a;(°)(0) = 1, y'-°HO) = or a;(°)(0) = 0, y'-^^O) 



(11) 



Equations (ITo| are subject to the initial conditions that 
1, while (11) is solved with the conditions x^^^(O) 



0. 



Solving the linear equations for the two sets of ini- 
tial conditions above and evaluating the values of the 
functions after one period T = 27r, one will obtain the 
following Floquet matrix 



J 



M 



^7V7r(-l + 3A/'2) 1- 



N 



(12) 



Note that the perturbation series (|9| is nonuniform. The 
approximation above breaks down at z ~ a~^/^, which 
limits the period T for a given drive amplitude a. 
Calculating the eigenvalues of M, we obtain that 



Ai,2 = 1 ± NnaVl - STV^ _ (9(^2 



(13) 



Because A > 1 for < TV < l/>/3, to the leading order 
we conclude that the continuation of the symmetric state 
is unstable in the interval and stable otherwise. The sta- 
bility is in agreement with Fig. [4] above. The instability 



for N < 1/ ^/3 implies that the symmetry breaking bifur- 
cation occurs earlier than for the undriven case, which 
we also observed numerically (not shown here). The an- 
alytical approximation (13) is compared to the numerics 



in Fig. [8[a), where rather perfect agreement is obtained. 
To illustrate our flnding, we show in Fig. [S] the time dy- 
namics of the initial condition A(0) = —Na, 0{0) = for 
N = 2c = 0.4 and N = 2c = 1.4 with a 0.1 corre- 
sponding to unstable and stable cases, respectively. The 
same method can be applied to analyse the existence and 
stability of the continuation of equilibrium solutions dis- 
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(c) 



FIG. 8: (a) The real and imaginary part of the critical Floquet multipliers shown in black and red solid curves respectively as 
a function of the solution norm N with a — 0.1. The dashed curve is our approximation ( |13| . (b-c) The time-dynamics of the 
power imbalance A(z) and the phase- difference 0{z) for the initial condition A(0) = —Na, 0(0) = and (b) N — 2c — 0.4 and 
(c) = 2c = 1.4, corresponding to an unstable and stable symmetric state, respectively. 



cussed in the previous sections for \a\ <C 1. However, in 
general the explicit expression of the solutions as well as 
the Floquet matrices will be very lengthy. 

VII. CONCLUSION 

Through numerically solving parametrically driven 
coupled nonlinear Schrodinger equations describing the 
dynamics of wavefields in an oscillating double-well po- 
tential, we have shown that such parametric drives may 
stabilize or destabilize the continuations of equilibrium 
time-independent states that are respectively unstable or 
stable in the undriven case. The analysis is performed 
by employing numerical continuations to find periodic 



orbits when varying a parameter and by calculating the 
corresponding Floquet multipliers of the states. Analyt- 
ical calculations that accompany the numerical results 
using an averaging method and perturbation expansion 
have been presented where good quantitative agreement 
is obtained. 
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